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pLn ' Abstract. Particle production sources at RHIC and LHC energies are investigated in 

QQ I pseudorapidity space. A nonequilibrium-statistical rclativistic diffusion model (RDM) 

with three sources is applied to the analysis of chargcd-hadron distributions in AuAu 
collisions at RHIC energies, in PbPb collisions at the current LHC energy of 2.76 
(~| ' TeV, in pPb at 5.02 TeV, and in pp. The size of the midrapidity source relative to 

the fragmentation sources in heavy-ion collisions is investigated as function of incident 
pi |. energy. At LHC energies, the midrapidity value is mostly determined by particle 

5^ ' production from gluon-gluon collisions. 
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, 1. Introduction 

Charged-hadron production in relativistic heavy-ion collisions has been investigated in 
great detail at the Relativistic Heavy-Ion Collider RHIC in AuAu collisions, and more 
O ■ recently at the Large Hadron Collider LHC in PbPb collisions. In particular, high- 
precision pseudorapidity distributions dNch/dr] of produced charged particles including 
their centrality dependence are now available in an energy range from ^/SNN = 0.019 
^ ■ to 2.76 TeV [H [2]- At RHIC energies these data include the fragmentation regions up 
to the values of the beam rapidities, whereas at the current LHC energy of 2.76 TeV 
corresponding to a beam rapidity of 7.99 preliminary but very precise ALICE data are 
available up to r/ ~ 5 [2J. 

Theoretical descriptions of the underlying partonic processes often focus on gluon- 
gluon production, such as in many approaches based on the color glass condensate 
(see |3] as an example). Based on this mechanism particle and antiparticle distributions 
would, however, be identical - which is not the case experimentally, as found for example 
in 77+ and tt~ distribution functions 

The relevance of the fragmentation sources from quark-gluon interactions has been 
investigated in a recent QCD-based study of net-baryon distributions (baryons minus 
antibaryons) . There the gluon-gluon source that is peaked at midrapidity cancels out 
such that only the fragmentation sources remain |5l |6], giving rise to two fragmentation 
peaks that are clearly seen in the data at high SPS and RHIC energies, and in the 
theoretical predictions at LHC energies. At low SPS energies the fragmentation peaks 
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overlap in rapidity space and hence, are not directly visible in the data, but can still be 
extracted quite reliably [7]. 

For produced particles (rather than net baryons), the effect of the fragmentation 
sources is less obvious, but clearly has to be considered. In this note I propose 
to investigate the relative importance of gluon-gluon vs. fragmentation sources as 
a function of cm. energy in collisions of heavy systems (AuAu, PbPb) using a 
phenomenological nonequilibrium-statistical model. This relativistic diffusion model 
(RDM) [8] has proven to be useful in the analysis of data and in predictions for 
asymmetric [9j and symmetric [10] systems, and - in spite of the small transverse size of 
the system - also in pp collisions [H] . Its three sources correspond to the gluon-gluon and 
fragmentation sources of the available microscopic theories. In direct comparisons with 
the available data the RDM can be used to infer the relative sizes of these underlying 
components as functions of the incident energy. 

In charged-hadron production at SPS and low RHIC energies up to i/J/wv — 20 
GeV, the gluon-gluon source centered at midrapidity is expected - and has turned 
out - to be unimportant [12], and the measured pseudorapidity distributions are well 
reproduced from the fragmentation sources only. At these relatively low energies, 
the fragmentation sources are peaked close to midrapidity and hence, are influenced 
considerably by the Jacobian transformation from rapidity to pseudorapidity space. 
At higher energies, the fragmentation peaks move apart, and the central gluon-gluon 
source emerges. Then the Jacobian increasingly affects only the central source. Also, its 
overall effect becomes smaller with rising energy since it depends on {{rn)/pTY. Still, 
a precise determination of the Jacobian is essential for the modeling of pseudorapidity 
distributions at LHC energies. The pronounced midrapidity dip that is seen in the 
recent ALICE PbPb charged-hadron data is due to the interplay of fragmentation and 
central sources, plus the effect of the Jacobian on the central source. 

A brief outline of the method used to determine the relative size and extent of the 
sources in r^— space is given in the next section. Results for heavy systems at RHIC 
and LHC energies are presented in section 3. The energy dependence of central and 
fragmentation sources is discussed in section 4. A brief outlook on pPb at 5.02 TeV and 
pp at 0.9 - 14 TeV is given in section 5, and the conclusions are drawn in section 6. 

2. Three sources model 

In the three-sources version of the relativistic diffusion model, rapidity distributions of 
produced particles are calculated from an incoherent superposition of the fragmentation 
sources Ri^2{y,t = Tint) with charged-particle content N^i^ (projectile-like), N^j^^ (target- 
like) and the midrapidity gluon-gluon source Rgg{y,t = Tint) with charged-particle 
content A^fif as 



dNch{y,t = Tint) 

dy 



'iMy, Tint) + Nl^R2{y, Tint) + N'^^Rgg{y, Tint) (1) 
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Figure 1. (Color online) The RDM pseudorapidity distribution function for 
charged liadrons in 0-5% central PbPb collisions at LHC energies of 2.76 TeV is 
shown in the upper panel of the figure, with the RDM parameters adjusted in a 
X^— minimization to the ALICE data (central value from [13], distribution from 
[2])- The fit takes the limiting fragmentation scaling hypothesis into account [lOj . 
The corresponding RDM parameters are given in table [1] In the lower part of the 
upper frame, calculated pseudorapidity distributions of produced charged particles 
from AuAu collisions (bottom) at ^^snn = 0.13 and 0.2 TeV for 0-6% central collisions 
with PHOBOS data [Tj are shown for comparison, see [M]. In the bottom frame, the 
underlying theoretical distributions are shown for 2.76 TeV PbPb. Only the shape of 
the midrapidity source is modified by the Jacobian. At LHC energies, the midrapidity 
value is mostly determined by particle production from gluon-gluon collisions. The 
upper curve is the RDM-prediction for 5.52 TeV. 





with the rapidity y = 0.5 ■ \n{{E + p)/{E — p)), and the interaction time r^nt (total 
integration time of the underlying partial differential equation). 

In the linear version of the RDM [8], the macroscopic distribution functions are 
solutions of the Fokker-Planck equation (A; = 1,2,3) 



d_ 

di 



Rk(,y,t) 



^ydy 



iveg - y) ■ Rk{y,t) 



+ D 



dy' 



■Rk{y,t)- 



(2) 
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The use of the Lorentz-invariant variable rapidity in the nonequihbrium-statistical 
Fokker-Planck framework has proven to be a useful approach in calculations and 
predictions of macroscopic distribution functions for produced particles. Integrating 
the equation with the initial conditions Ri^2{y, t = 0) = S{y ± Umax), the absolute value 
of the beam rapidities ymax, and R3=gg{y,t = 0) = 6{y — yeq) yields the exact solution. 
The mean values are derived analytically from the moments equations as 

< yi,2{t) >= Vegl^ - exp{-t/Ty)] =F ymax {-t/Ty) (3) 

for the sources (1) and (2) with the absolute value of the beam rapidity ymax and the 
rapidity relaxation time Ty. 

The local equilibrium value yeq is equal to zero only for symmetric systems. For 
asymmetric systems such as pPb, the midrapidity source is moving [15], and the 
superposition of the sources is even more sensitive to the values of the model parameters 
than in the symmetric case. From energy-momentum conservation the centrality- 
dependent equilibrium value is obtained as [HI [17] 

(6) = -0 5 ■ In ("^n^)) exp(y^a^) + (m^(&)) exp{-ymax) 

exp ( ymax ) 

with the beam rapidities ybeam = TVmax, the average transverse masses {mi2{b)) = 
\Jrni 2{b) + (pt)'^, and participant masses mi 2(6) of the p- and Pb-like participants in 
pPh collisions that depend on the impact parameter b. The minus sign refers to cases 
where > such as in the ALICE pPh experiment of 2012 where the Pb beam 
defined the positive rapidity. The sign of the equilibrium value changes when the beams 
are interchanged, as is planned in the 2013 pPb experiments at the LHC. 

For sufficiently large beam rapidities ymax such at LHC energies, the equilibrium 
value can be approximated as 

,.,(6). 0,5.. „ nil . (5) 

The corresponding numbers of participants can be obtained from the geometrical 
overlap, or from Glauber calculations. For yteam — 8.6 as in 5.02 TeV pPh an estimate 
in — 5% central collisions is yeq — 0.6, and smaller values for more peripheral collisions. 
The time evolution in the RDM causes a drift of the distribution functions Ri^2 towards 

Veq- 

Whether the mean values of Ri and R2 actually attain yeq depends on the centrality- 
dependent interaction time Tmt (the time the system interacts strongly, corresponding 
to the integration time of ([2])), and its ratio to the rapidity relaxation time Ty. Typical 
interaction times at LHC energies from dynamical models in central PbPb collisions are 
6-8 fm/c, which is too short for the fragmentation sources to reach equilibrium, such 
that their mean values < yi^2 > remain between beam and equilibrium values. 

This does not apply, however, to Rgg which already emerges near equilibrium at 
the parton formation time - here, at t = because of the (idealized) 5— function initial 
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condition - and spreads in time due to strong diffusive interactions with other particles, 
without any shift in the mean value for a given centrality class. The variances are 

crl{t) = D';Ty[l~eM-'^t/Ty)], (6) 

so they reach equilibrium faster than the mean values. Here the diffusion coefficients in 
rapidity space are D^, and presently I assume equal values for the three sources. The 
corresponding FWHM- values Ti^2,gg are listed in table H] for AuAu and PbPb together 
with the mean values. Both are determined from x^— minimizations with respect to 
the very precise PHOBOS and preliminary ALICE data [HE]- 

Since the theoretical model is formulated in rapidity space, one has to transform 
the calculated distribution functions to pseudorapidity space, rj = — ln[tan(6'/2)], in or- 
der to be able to compare with the available data, and perform x^— minimizations. The 
well-known Jacobian transformation 



dN dNdy ^, , .dN 

-r = iri~ = •^^'^^^'P^'IT' (7) 
drj dy drj dy 

J{r],m/pT) = cosh(?7)-[l + {m/pTY + sinh^(?7)]^^/^ (8) 

depends on the squared ratio of the mass and the transverse momentum of the produced 
particles. Hence, its effect increases with the mass of the particles, and it is most 
pronounced at small transverse momenta. For reliable results one has to consider the 
full pr— distribution, however: It is not sufficient to consider only the mean transverse 
momentum {pt)- In [ID] we have discussed in some detail how this can be done for 
known pj-— distributions of identified ti~,K~^ and antiprotons. We use the pion mass 
m^, and then calculate an effective mean transverse momentum < pjf'^ > such that the 
experimentally determined Jacobian Jy=Q of the total charged-hadron distribution at 
rapidity zero is exactly reproduced. This yields for a given centrality class [10] 

{p'Jf) = m.J,=o/v/l - JU ■ (9) 
These effective transverse momenta are smaller than the mean transverse momenta 
determined from the pr— distributions, and the corresponding effect of the Jacobian is 
therefore larger than that estimated with {p^) taken from the transverse momentum 
distributions for each particle species. At high RHIC and LHC energies the effect of 
the Jacobian transformation remains, however, essentially confined to the midrapidity 
source. 

The Jacobians can now be calculated for each centrality class, pseudorapidity 
distributions of produced charged hadrons are obtained in the three-sources model from 
([T]), the parameters are optimized with respect to the available data, and conclusions 
regarding the relative sizes of the sources become possible. 

However, LHC data are still missing in the fragmentation region. We have therefore 
proposed in [TO] to use the well-known limiting fragmentation scaling hypothesis [18] 
as an additional constraint: At sufficiently high energy, particle production in the 
fragmentation region becomes almost independent of the collision energy. Hence we use 
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Figure 2. (Color online) Pseudorapidity distributions for produced charged hadrons 
in 2.76 TeV PbPb collisions as functions of centrality, from bottom to top: 20-30%, 10- 
20%, 5-10%, 0-5%. Calculated RDM distributions (solid curves) have been optimized 
in x^-fits with respect to the preliminary ALICE data from [2] , and using the limiting 
fragmentation scaling hypothesis in the region of large rapidities where no data are 
available. The centrality-dependent parameter values are as in [TU] . 



0.2 TeV AuAu results at RHIC - where data in the fragmentation region are available 
- to supplement the LHC 2.76 TeV PbPb data in analogous centrality classes at large 
values of pseudorapidity, shifting the latter by Ay = ybiim~ybean? — 7.99 — 5.36 = 2.63. 
The resulting RDM-parameters have physically reasonable dependencies on the cm. 
energy and centrality, and their extrapolations - in particular, to higher energies such 
as the LHC design energy of 5.52 TeV PbPb - can readily be used for predictions. In this 
work the main emphasis is on the determination of the relative source sizes as function 
of the incident energy. For a reliable identification of the percentage of particles in the 
three sources the accurate modeling of the fragmentation region is essential, and limiting 
fragmentation scaling is a useful tool for PbPb collisions at LHC energies since accurate 
data at RHIC energies are available for AuAu at large values of pseudorapidity. 

For pPb collisions, however, no low-energy data in the fragmentation regions exist, 
and limiting fragmentation scaling can not be used at LHC energies. Also, the presently 
available ALICE data at 5.02 TeV are limited to the range -2 < r/ < 2 pjj. It is 
therefore advisable to wait for the forthcoming 2013 LHC pPb data in the pseudorapidity 
region 2 < < 5 before more definite conclusions regarding the relative source sizes can 
be drawn. 

For pp at LHC energies there is presently also a lack of data in the fragmentation 
region. The pp data at 0.9, 2.36 and 7 TeV [201 HI] have a pseudorapidity coverage 
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\ri\ < 2.25. There exist TOTEM inelastic results at 7 TeV for pseudorapidities 5.4- 
6.4 [22], but these are still far from the beam rapidity, and their normalization is not 
consistent with the non-single diffraction (NSD) CMS results at midrapidity. As a 
consequence, the determination of the relative source sizes as function of the incident 
energy is restricted in this work to heavy-ion collisions, although I shall also present a 
pPh analysis at 5.02 TeV, and pp RDM results from the earlier work [H] in comparison 
with the recent TOTEM data. 

3. Results 

The result of the three-sources RDM calculation for the pseudorapidity distribution of 
produced charged hadrons 2.76 TeV PbPb is shown in figure [1] together with recent 
preliminary ALICE data [2j for — 5% centrality in a optimization. Parameters are 
given in table [H The published midrapidity ALICE data point [13] is slightly below, 
but within the error bars compatible with the more recent data. 

The three-sources fit uses the limiting fragmentation scaling assumption based on 
the 0.2 TeV AuAu central PHOBOS data [T] from RHIC, which are also shown in figure 
[D together with the 0.13 TeV AuAu data, and the corresponding RDM results [H], with 
parameters given in table [H Results for central AuAu collisions at RHIC energies of 
19.4 GeV and 62.4 GeV that are included in table [1] are taken from earlier work [TO] [T^. 
A prediction for the LHC design energy of 5.52 Te PbPb is also shown in figure [B 

The relative size of the three sources in central 2.76 TeV PbPb is shown in the 
lower frame of figure [H At this LHC energy, the midrapidity source already contains 
the largest fraction of produced charged hadrons. Its shape is significantly deformed 
by the Jacobian transformation from rapidity to pseudorapidity space, whereas the 
fragmentation sources are not much influenced by the transformation. 

In the full distribution that arises from the incoherent superposition of the three 
sources, it is evident that the midrapidity dip is more pronounced at LHC energies as 
compared to RHIC energies, although the effect of the Jacobian tends to be smaller at 
the higher incident energy. This clearly indicates that there has to be a physical origin 
of the midrapidity dip in addition to the effect of the Jacobian. 

The hypothesis promoted in this work is that the interplay of the three sources 
provides the observed effect. In 2.76 TeV PbPb collisions, the fragmentation sources 
are peaked at large values (< ?/i 2 >= 3.34) of rapidity - whereas at 0.2 TeV RHIC 
energy, the center is at < yi 2 >= 2.4. Consequently, the midrapidity yield at LHC 
energies is essentially due to the central source, with only a small contribution from 
the fragmentation sources. Although the relative particle content in the central source 
is larger at LHC energies than at RHIC, this produces the observed midrapidity dip, 
together with the effect of the Jacobian on the central source. 

The centrality dependence of charged hadron production at LHC energies as 
displayed in figure [2] in comparison with preliminary ALICE PbPb data for 0-5%, 5- 
10%, 10-20% and 20-30% shows that the total number of produced charged hadrons 
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Figure 3. (Color online) Number of produced charged hadrons as function of the cm. 
energy ^snn from RDM-fits of the available data for central heavy-ion collisions at 
0.019, 0.062, 0.13, 0.2 TeV (RHIC, AuAu), and 2.76 TeV (LHC, PbPb). Circles are 
the total numbers, following a power law oc s^^. Squares are hadrons produced from 
the midrapidity source, with a power law cx s^^, and triangles are particles from the 
fragmentation sources cx log(sjvA'/'So)- The gluon-gluon source (dashed) becomes the 
main source of particle production between RHIC and LHC energies. 

rises with increasing centrality. As displayed in figure 2 of [10] tliat was based on 
earlier preliminary ALICE data, the rise is almost linear with increasing number of 
participants. The percentage of particles in the midrapidity source depends weakly on 
centrality, falling from 56 % in central collisions to 51 % at 20-30 % centrality. 

4. Energy dependence of the hadron production sources 

There are now sufficiently precise data on charged-hadron production at RHIC [1] and 
LHC p] energies available in order to investigate the relative size of the three particle 
production sources as function of energy in heavy-ion collisions (AuAu at RHIC, PbPb 
at LHC). Based on optimizations of the analytical three-sources RDM solutions with 
respect to these data I have displayed the energy dependence of the sources in figure |3l 
with parameters as shown in table [H 

According to these results, the total charged-hadron production (circles) follows a 
power law oc s^^. The hadrons produced from the central source (squares) have an even 
stronger dependence on initial energy according to oc s^'^, whereas particles produced 
in the fragmentation sources have a weaker dependence oc log(sArAr/so). 

The strong rise of the particle production yield from the central (gluon-gluon 
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Table 1. Three-sources RDM-parameters Tint/ry, ri^2,rgg, and Ngg. N^^"^ is the 
total charged-particle number in the fragmentation sources, Ngg the number of charged 
particles produced in the central source. Results for < 1/1^2 > are calculated from 
Vbeam and Tint/Ty. Valucs are shown for 0-5% PbPb at LHC energies of 2.76 and 
5.52 TeV in the lower two lines, with results at 2.76 TeV from a x^-minimization 
with respect to the preliminary ALICE data [2 , and using limited fragmentation as 
constraint. Corresponding parameters for 0-6% AuAu at RHIC energies are given for 
comparison in the upper four lines based on PHOBOS results [1]. Parameters at 5.52 
TeV denoted by * are extrapolated. Experimental midrapidity values (last column) 
are from PHOBOS [1] for |??| < 1, 0-6% at RHIC energies and from ALICE [13] for 
|?7| < 0.5, 0-5% at 2.76 TeV. 
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induced) source is evidently due to the increasing gluon content of the system at high 
relativistic energies. In particular, the total particle production rate from the central 
source becomes larger than that from the two fragmentation sources at an incident 
energy between the highest RHIC energy (0.2 TeV), and the LHC regime. In view 
of the lack of data in this intermediate regime, the precise crossing point is, however, 
difficult to determine. 

5. Proton-induced collisions 

I had previously applied the relativistic diffusion model [8] also to pp collisions at RHIC 
and LHC energies [TT] . Here transport phenomena are not expected to be fully developed 
because of the small transverse size of the system, but still, the RDM yields reasonable 
results for pseudorapidity distributions of produced charged hadrons in the 77— range 
where data are available, and provides predictions for large values of 77. 

The TOTEM collaboration has reported first experimental pp results at 7 TeV 
and large |r/| [22]. Measured values for dN/drj of charged hadrons range from 
3.84±0.01(stat)±0.37(syst) at \r]\ = 5.375 to 2.38± 0.01(stat)±0.21(syst) at \ri\ = 6.375. 
The data account for about 95% of the total inelastic cross section. A corresponding 
RDM prediction taken from [TT] is compared with these new data in figure HI It has 
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Figure 4. Pseudorapidity distributions of produced charged hadrons in pp collisions 
(NSD) at LHC cm. energies of 0.9, 2.36, 7 and 14 TeV (bottom to top) as calculated 
in the three-sources approach [ll] and fitted to CMS NSD data (20l E]. At 0.9 TeV 
UA5 NSD data are also shown |23], triangles. Data points in the region 77 — 5.375 
to 6.375 are results from the TOTEM collaboration [22] at 7 TeV (95% inelastic; the 
prediction was normalized to NSD). 

the correct slope, whereas MC predictions give a different slope The absolute 

magnitude of the RDM result is, however, slightly too low since it was normalized to 7 
TeV CMS NSD pp data at midrapidity, which are below inelastic results. 

In order to draw definite conclusions regarding the particle content of the three 
sources in pp collisions as function of ^/s one would need data at various incident energies 
over a larger range in pseudorapidity than what is presently available. The same is true 
for charged-hadron production in asymmetric proton-induced collisions. Here ALICE 
data from the 2012 testrun are available for pseudorapidity distributions of produced 
charged hadrons in pPh collisions at y/SpN = 5.02 TeV (corresponding to a proton beam 
momentum of 4 TeV/c) in the range —2 < t] < 2 [19j . 

The midrapidity yield is almost two orders of magnitude smaller than in 0-5% 
central PbPb collisions at 2.76 TeV, and the shape of the distribution function is very 
sensitive to the details of the underlying partial distributions, see figure [51 The central 
source is seen to be significantly modified by the Jacobian, whereas the asymmetric 
fragmentation sources are almost gaussian-shaped. 

The interplay of the fragmentation sources with the moving gluon-gluon source 
that is centered at small positive 77 values generates the characteristic shape of the 
charged-hadron distribution. The proton-like side of the distribution is seen to be 
considerably steeper than the Pb-like side. It is expected that the 2013 pPh runs with 
interchanging beams will enable the LHC heavy-ion experiments to actually measure this 
predicted difference. A corresponding effect had been observed in dAu at 0.2 TeV [HE]. 
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Figure 5. (Color online) The predicted RDM pseudorapidity distribution function 
for charged hadrons in minimum bias pPh collisions at LHC cm. energy of 5.02 TeV 
shown here is adjusted in the mid-rapidity region to the preliminary ALICE data |19j 
(systematic error bars only). The underlying distributions in the three-sources RDM 
are also shown, with the dashed curve arising from gluon-gluon collisions, the dash- 
dotted curve from valence quark-gluon events in the Pb-like region, and the dotted 
curve correspondingly in the proton-like region. 



There it was also shown that the difference in the slopes becomes more pronounced 
with increasing centrality, and a similar behaviour is expected at LHC energies once 
centrality-dependent data become available. 

For a reliable determination of the particle content in the three sources one needs 
data at larger values of pseudorapidity, which are not yet available. Moreover, for pPh 
only a single incident energy is available, so that an investigation of the relative particle 
content in the three sources as function of y^Sp^ similar to the heavy-ion case is presently 
not feasible. 

6. Conclusions 

The particle content of fragmentation and midrapidity (gluon-gluon) sources for 
charged-hadron production in heavy-ion collisions at high relativistic energies has been 
determined as function of cm. energy in a phenomenological approach. Due to the 
availability of preliminary high precision PbPb ALICE data from LHC [2], earlier 
AuAu PHOBOS data [1], and with the limiting fragmentation scaling hypothesis in 
the region of large pseudorapidities where LHC data are still missing, a determination 
of the particle content in a three-sources non-equilibrium statistical model has become 
possible. 
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Whereas the fragmentation sources are found to depend on energy oc log(sArAr/so), 
the central source has a much stronger energy dependence oc s^"^, and the total number 
of produced charged hadrons - which arises from an incoherent superposition of the 
three sources - behaves like oc s^^. As a consequence, particle production from the 
gluon-gluon source becomes more important than that from the fragmentation sources 
in the energy range between the maximum RHIC energy of 0.2 TeV, and the current 
LHC energy of 2.76 TeV. 

The same approach has also been applied to charged-hadron production in pp 
collisions at LHC energies. This yields predictions for the pseudorapidity distributions 
of produced charged hadrons in the region of large rj which are in reasonable agreement 
with recent inelastic data of the TOTEM collaboration [22] for 7 TeV pp, although there 
is a normalization problem when comparing with the midrapidity NSD results. Due to 
the lack of large-r] data at other LHC energies, a precise determination of the particle 
content of the sources as function of cm. energy is, however, presently not yet feasible 
in proton-proton collisions, although there exist very precise PHOBOS pp data at RHIC 
energies of 0.2 and 0.41 TeV [1] that extend into the fragmentation region. 

In central pPh collisions at the LHC cm. energy of 5.02 TeV, ALICE data [19] have 
been used to compare with the analytical solutions of the relativistic diffusion model 
(RDM). The shape of the pseudorapidity distribution is found to be very sensitive to the 
interplay of the p— and Pb-like fragmentation sources, and the moving central source 
that is significantly modified by the Jacobian. The RDM calculation exhibits a steeper 
slope on the proton-like side, as compared to the Pb-like side. It is expected that the 
forthcoming LHC pPh experiments confirm this behaviour. 
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